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1,  INTRODUCTION 


The  prediction  of  finned  kinetic  energy  (KE)  projectile  base  flow  is  a  challenging  problem 
in  the  area  of  applied  computational  fluid  dynamics  (CFD).  Finned  projectile  base  flow, 
even  for  a  non-spinning  configuration  at  0°  angle  of  attack,  is  a  complex  3-D  flow  field  whose 
modeling  requires  the  most  powerful  supercomputing  resources  available.  Accurate  modeling 
of  the  base  flow  will  provide  the  base  drag  coefficient,  one  of  the  remaining  aerodynamics 
coefficients  yet  to  be  adequately  predicted  for  KE  projectiles  using  CFD.  A  CFD  capability 
for  KE  projectile  base  flow  will  provide  an  understanding  of  the  flow  structure,  possibly 
leading  to  the  design  of  effective  drag-reduction  techniques  and  lethality  enhancements  for 
finned  projectiles. 

This  report  documents  a  CFD  study  of  the  base  flow  of  a  105-mm  M735  finned  KE 
projectile.  Two  Navier-Stokes  finite  difference  techniques  are  used  in  conjunction  to  generate 
flow  field  solutions  at  Mach  numbers  3,  4,  and  5.  The  major  objectives  of  the  study  are  to 
(1)  numerically  simulate  the  projectile  base  flow,  (2)  compare  the  total  computed  drag  of 
the  projectile  with  range-acquired  data,  (3)  gain  insight  into  the  flow  field  structure,  and  (4) 
characterize  the  computational  requirements  of  the  problem. 

An  informative  overview  of  projectile  CFD  work  at  the  U.S.  Army  Research  Laboratory^ 
(ARL)  is  provided  by  Sturek,  Nietubicz,  Sahu,  and  Weinacht  (1992).  That  report  traces  the 
development  and  application  of  CFD  methods  for  computing  KE  projectile  aerodynamics 
(excluding  the  base  flow),  as  well  as  axisymmetric  and  3-D  base  flows  for  bodies  of  revolution. 
The  present  study  represents  the  first  in-house  attempt  to  model  the  base  flow  of  a  finned 
projectile  using  CFD. 


2.  CONFIGURATIONS  AND  FLIGHT  CONDITIONS 

The  configuration  of  interest  is  the  M735  projectile,  which  is  a  fin-stabilized,  long-rod  KE 
projectile  fired  from  a  105-mm  gun  tube.  The  computational  model  of  the  M735  projectile, 
illustrated  in  Figure  1,  is  a  simplified  version  of  the  actual  model.  The  model  possesses 
a  conical  nose  section,  followed  by  a  cylindrical  section  of  constant  reference  diameter,  d, 
equal  to  35.2  mm.  The  total  model  length  is  13.94  calibers  (1  caliber  =  1  cal  =  1  reference 
diameter).  The  blunt  nosetip  is  replaced  by  a  conical  nosetip,  and  the  sabot  grooves  are 
replaced  by  a  smooth  surface. 


*  Formerly  the  U.S.  Army  Bsdlistic  Research  Laboratory  (BRL) 
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The  six  fins  are  equally  distributed  around  the  body  and  aligned  with  the  projectile 
axis.  The  leading  edges  are  cylindrically  blunted  and  swept  71°.  The  trailing  edges  are 
truncated,  perpendicular  to  the  projectile  axis,  and  aligned  with  the  projectile  base  (i.e.,  no 
fin  overhang).  The  fins  are  tapered  from  root  to  tip,  as  shown  in  Figure  1.  The  leading  and 
trailing  edge  chamfers,  which  control  the  spin  of  the  actual  projectile,  are  omitted  from  the 
computational  model.  The  juncture  of  the  cylinder  and  the  fin  leading  edge  is  modified  with 
a  small  rounded  fillet,  which  typically  extends  about  0.3  to  0.5  cal  upstream  and  downstream 
from  the  juncture.  The  base  of  the  computational  model  is  assumed  to  be  perfectly  flat,  in 
which  protuberances  such  as  the  tracer  cavity  and  the  core  rod/finned  afterbody  interface 
are  ignored. 

Computational  results  are  presented  for  free  stream  Mach  numbers  of  3,  4,  and  5.  The 
angle  of  attack  is  fixed  at  0°  and  the  spin  rate  is  fixed  at  0  rpm.  The  free  stream  conditions 
are  taken  to  be  sea-level  atmospheric  values.  The  projectile  surface  temperature  is  specified 
as  294  K  or  adiabatic,  and  is  discussed  in  detail  later.  Fully  turbulent  boundary  layer  flow 
is  assumed  over  the  projectile  forebody  and  fins.  Laminar  flow  is  assumed  in  the  base  region 
for  this  initial  computational  effort. 

The  roll  angle,  of  a  point  in  the  flow  or  on  the  body  is  defined  as  the  circumferential 
angle  as  measured  relative  to  one  reference  fin,  as  illustrated  in  Figure  2.  Since  the  angle  of 
attack  IS  0°  and  the  configuration  has  six  fins,  only  roll  angles  between  0°  and  30°  are  nec¬ 
essary  for  clear  reference  because  of  circumferential  symmetry  considerations.  For  practical 
purposes,  the  convention  is  henceforth  adopted  that  =  0°  refers  to  any  of  the  six  circum¬ 
ferential  planes  {<t>  =  60°  or  ^  =  180°,  for  example)  that  are  coplanar  with  a  fin  mid-plane, 
and  4>  =  30°  refers  to  any  of  the  six  circumferential  planes  that  are  located  midway  between 

two  adjacent  fins.  Mirror  symmetry  in  the  geometry  and  flow  field  is  assumed  to  exist  at 
roll  angles  (f)  =  0°  and  =  60°. 


3.  FLOW  MODEL  FOR  PROJECTILE  FOREBODY  AND  FINS 

The  complete  flow  field  surrounding  the  projectile  is  computed  using  two  separate  3-D 
Navier-Stokes  numerical  techniques  in  conjunction,  as  shown  in  Figure  3.  The  first  is  a  parab¬ 
olized  Navier-Stokes  (PNS)  technique,  which  computes  the  flow  over  the  projectile  forebody 
and  fins.  The  second  is  an  unsteady  Navier-Stokes  (UNS)  technique,  which  computes  the 
projectile  base  flow.  The  PNS  technique  is  briefly  described  next;  the  UNS  technique,  the 
application  of  which  is  the  focus  of  this  research  effort,  is  discussed  in  some  detail  thereafter. 
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3.1  Parabolized  Navier-Stokes  CFD  Technique.  As  shown  in  Figure  3,  the  PNS 
technique  is  used  to  provide  the  upstream  inflow  boundary  condition  for  the  base  flow  compu¬ 
tation.  The  PNS  technique  also  provides  the  forebody  drag  contribution  (from  the  forebody 
and  fins),  which  is  added  to  the  base  drag  contribution  (from  the  cylinder  base  and 

fin  trailing  edge),  to  obtain  the  total  computed  zero-yaw  drag  coefficient,  Cdqi  that  is 
compared  to  range-acquired  values.  The  values  are  tabulated  in  the  results  section. 

The  PNS  technique  has  been  adapted  and  extensively  used  within  the  Weapons  Technol¬ 
ogy  Directorate  of  ARL.  The  original  technique  was  developed  by  Schiff  and  Steger  (1979), 
and  the  first  in-house  application  to  U.S.  Army  projectiles  was  done  by  Sturek  and  Schiff 
(1981).  Applications  of  the  technique  to  KE  and  other  projectiles  have  been  documented 
in  reports  by  Rai,  Chaussee  &  Rizk  (1983);  Weinacht,  Guidos,  et  al  (1985,  1986,  1993); 
Weinacht,  et  al  (1988,  1990,  1991);  and  Guidos  (1994).  PNS  techniques  are  space-marching 
(as  opposed  to  UNS  techniques,  which  are  time-marching),  in  which  one  numerical  integra¬ 
tion  sweep  is  made  from  the  nosetip  of  the  projectile  to  the  beise  to  obtain  a  single  steady 
state  solution.  Because  of  this  space-marching  formulation,  PNS  techniques  are  orders  of 
magnitude  more  computationally  efficient  than  UNS  techniques  for  computing  the  attached 
regions  of  supersonic  flow  of  interest  here.  The  PNS  solutions  for  this  configuration  were 
generated  using  about  1  to  2  hours  of  central  processing  unit  (CPU)  time  per  solution  on  a 
Cray  X-MP  computer. 

This  PNS  technique  is  fully  implicit,  approximately  factored  (Beam  &  Warming,  1978), 
and  uses  second-order  central  differencing.  In  this  formulation,  only  the  thin-layer  viscous 
terms  in  the  outward  coordinate  direction  are  retained.  The  Baldwin-Lomax  (1978)  two-layer 
algebraic  turbulence  model  is  used  in  all  cases  presented  here.  In  addition,  a  shock-fitted 
outer  boundary  condition  reported  by  Rai  and  Chaussee  (1983)  is  used.  The  code  is  used  in 
step-back  mode  to  generate  a  pointed  conical  nosetip  solution,  providing  the  initial  conditions 
for  downstream  marching,  as  discussed  by  Sturek  and  Schiff  (1981). 

Typical  of  long-rod  KE  projectiles,  the  actual  projectile  possesses  sabot  grooves  over  the 
cylindrical  forebody  section.  These  grooves  act  to  transfer  the  accelerating  force  from  the 
sabot  (which  also  acts  as  a  borerider)  to  the  projectile.  The  correlation  of  Mikhail  (1989)  was 
used  here  to  estimate  the  effect  of  groove  drag  for  the  M735  configuration.  The  correlation 
shows  that  the  sabot  grooves  account  for,  at  most,  1.6%  of  the  total  drag  for  Mach  numbers 
between  3  and  5.  With  this  result  in  hand,  all  computational  drag  results  presented  in  this 
report  exclude  any  contributions  from  the  sabot  grooves.  As  already  mentioned,  the  grooves 
are  replaced  by  a  smooth  surface  in  the  computational  model. 
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Two  sets  of  PNS  computations  were  generated.  The  first  set  specified  the  projectile 
surface  temperature  to  be  294  K.  This  condition  is  representative  of  the  temperature  of  the 
projectile  in  early  flight,  used  for  obtaining  the  viscous  contribution  to  the  forebody  drag. 
The  second  set  of  PNS  computations  specified  the  projectile  surface  temperature  condition 
to  be  adiabatic,  used  as  the  inflow  condition  for  the  base  flow  computations. 

4.  FLOW  MODEL  FOR  PROJECTILE  BASE  REGION 

4.1  Unsteady  Navier-Stokes  CFD  Technique.  The  UNS  technique  is  a  time¬ 
marching  scheme,  in  which  an  initial  assumed  flow  field  is  integrated  through  time  toward 
a  steady  state  solution.  The  technique  used  here  has  previously  been  described  and  applied 
to  hypersonic  vehicles  (Edwards,  Chaussee,  Lawrence,  &  Rizk,  1987;  Edwards,  1988)  and 
referred  to  as  the  UWIN  code.  The  version  of  the  code  used  here  is  written  in  a  two-zone 
grid  formulation,  as  illustrated  in  Figure  3. 

4.2  Governing  Equations.  The  governing  equations  are  a  nondimensionalized  form 
of  the  time-dependent,  Reynolds- averaged,  thin-layer  Navier-Stokes  equations  in  3-D  gener¬ 
alized  coordinates  and  strong  conservation  law  form,  given  as  (Rai  1989) 


Qt  +  +  Fj,  +  (j(;  =  Rj  ^5^  (1) 

The  generalized  curvilinear  coordinates  are 
T  =  t  is  time 

C  is  the  longitudinal  (or  ^streamwise”  or  *^near- axial”)  coordinate 

T]  =  r]{x,y,z^t)  is  the  near-circumferential  coordinate 
C  =  is  the  outward  coordinate 


The  vector  of  dependent  variables,  is 

Q  =  JQ  =  [p^  pu,  pv,  pw,  ef  (2) 

in  which  the  density  is  p;  the  velocity  components  in  the  a:,  y,  and  z  directions  are  u,  u,  and 
respectively;  and  the  total  energy  per  unit  volume  is 
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The  vectors  E,  F,  and  G  are  the  transformed  inviscid  flux  vectors: 


-  T-i 


E  =  J 


pU 

puu + ixp 
pvU  rf-  tyP 
pwU  + 

{e  +  p)U 


F  =  J- 


pV 

puV  +  TjxP 
pvV  +  r}yP 
pwV  +  •qzp 
{e-\-p)V 


G  =  J-i 


pW 

puW  +  CrP 
pvW  +  CvP 
pwW  +  Gp 
{e  +  p)W 


in  which  U,  V,  and  W  are  the  velocities  in  the  transformed  coordinate  directions,  i.e., 

U  =  +  + 

V  =  T]xU  +  T]yV  + 

+  CyV  +  CzW  (4) 

Thin-layer  viscous  terms  are  retained  only  in  the  (  direction  and  are  contained  in  the 
vector  S,  given  as 

r  0 

TniU(  -f  m2Cx 

S  =  J~^  miV(  -f  m2Cy  (5) 

miu;^  -1-  m2Cz 


in  which 


mi  =  piC  +  Cy  +  C) 
m2  =  (p/3)(CxU<  + 

"•3  =  (ci+c+a^^(^)+|(«=+v=+^^)(i 

-t-m2(Cx«  +  CyV  +  Gw) 


The  following  nondimensionalization  is  used  to  define  the  dependent  variables: 


Plo  "  O^o/VT  OSo/v/7  ^  Plo^loVl 

in  which  the  superscript  “  *  ”  denotes  a  dimensional  quantity.  Perfect  gas  behavior  is 
assumed;  therefore,  the  pressure  is  defined  from  the  dependent  variables  as: 

P=  (7- 1)[£- +  +  (8) 

The  formulae  for  the  metrics  of  the  coordinate  transformation  (^t,  Gi  and  the 

Jacobian,  J ,  are  given  by  Pulliam  and  Steger  (1978).  In  regions  where  turbulence  is  modeled, 
the  Baldwin-Lomax  (1978)  model  is  used. 
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4.3  Numerical  Discretization  and  Solution  Algorithm.  The  governing  flow  equa¬ 
tions  are  discretized  using  an  approximately  factored,  implicit,  flux-split,  upwind,  total  vari¬ 
ation  diminishing  (TVD),  finite  difference  approach.  The  scheme,  used  here  to  generate 
steady  state  solutions,  is  temporally  first  order  accurate  and  spatially  second  order  accurate. 
Applying  the  scheme  to  the  governing  equations  yields  (omitting  most  unmodified  indices) 
the  following  discretization  (Rai  &  Chakravarthy  1986): 


’  + 


X 


At 


/  +  ^(v<C+  +  A<c--ij;'i;JV) 


X  {Q"*'  -  (}”) 


=  —At 


—  P^  P^  Pn 


At/ 


AC 


■^2AC  A£^''’((5j-i, Q|-, C)+i/2) 

+  AE-(g,,g,+i,c._v2)  -  AE-(g.+x,g,+2,c.+i/2)]" 

^2 At/  ^T-1/2)  ~  AF‘^((5j-i,  ^j,T/j+i/2) 


■h  at*  (Qj,  Qj+i,  T/j_i/2)  Ai^  Qi+2)  ^?j+i/2)] 

^2AC  ~  A (?■*■( 1,  Q;t,  Cfc-n/2) 


-h  AG  {Qk-iQk-^-iXk-i/i)  —  ((5fc+i,<?;b+2,Cfc-i-i/2)]  +  At 


On  _  On 

^fc-H/2  ^fc-1/2 

R^AC 


(9) 


in  which  terms  like  A^,  V,,,  and  S(^  are  forward,  backward,  and  central  difference  operators, 
respectively,  in  the  appropriate  coordinate  direction.  The  grid  point  indices  i,j,  and  k 
correspond  to  the  C,  ?/,  and  C  directions,  respectively.  The  superscript  n  represents  the 
current  time  level,  and  7  is  the  5x5  identity  matrix. 


The  implicit  (left-hand)  side  of  this  equation  has  been  approximately  factored  (Beam  & 
Warming  1978).  The  matrices  and  G±  are  the  inviscid  flux-split  Jacobian  matri¬ 

ces  which  relate  the  changes  in  the  inviscid  fluxes  with  respect  to  the  dependent  variables 
(Pulliam  and  Steger  1978),  evaluated  as 

=  d{E^)/dQ  =  d{F^)/dQ  =  d{G^)/dQ  (10) 

The  matrix  N  is  the  viscous  Jacobian  matrix  corresponding  to  the  C  direction,  whose  terms 
are  found  using  the  procedure  outlined  by  Steger  (1978). 
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The  explicit  (right-hand)  side  of  Eq.  (9)  includes  the  vectors  F,  and  G,  which  are 
first  order,  inviscid,  numerical  fluxes;  and  the  vectors  AF^,  and  AG^,  which  are 

flux  difl'erences  contributing  second  order  spatial  accuracy.  These  vectors,  key  elements  of 
the  numerical  scheme,  are  constructed  using  the  upwind  method  reported  by  Chakravarthy 
and  Osher  (1983).  The  flux  differences  are  defined  from  finite  differences  evaluated  along 
three  subpaths  in  a  state  space  defined  by  the  Riemann  invariants  that  correspond  to  the 
eigenvalues  of  the  Jacobian  matrices.  The  flux  differences  are  then  modified  using  a  flux- 
limiting  procedure  (Rai  and  Chakravarthy  1986),  and  denoted  as  AE'^,  AF'^,  and  AG^. 
Separately,  the  terms  in  the  viscosity  vector  S  are  numerically  differentiated  using  standard 
second-order  central  differences,  cis  given  by  Rai  (1989). 

Equation  (9)  is  a  linearized  form  of  the  nonlinear  governing  equations  presented  in  Eq. 
(1).  The  nonlinear  form  of  the  finite  difference  equations  is  solved  by  applying  a  Newton- 
iterative  technique  (Rai  and  Chakravarthy  1986).  The  algorithm  can  be  written  as 

I  +  X  !  +  ” 

X  /  +  -  R:'SiN)  X  (O'*'  -  Q’) 

A?  + - - + - AC - 

+  2^  [A.E^(Qi-2»Qi-ij 6-1/2)  —  AE^{Qi-i,Qi,(i+i/2) 

AE  ~  AE  (Q«+i,  Q«+2)  6+1/2)] 

"^2 At/  [A-^  {Qj-2iQj-ii‘0j-i/2)  ~  AF 

■b  A-F  (Qj,Qj+i,T}j_if2)  ~  AF  (Qj+1,  Qj+2,  T/j+.i/2)] 

^2AC  {Qk-2,Qk-\,Ck-l/2)  —  AG^ {Qk-l,Qk,Ck+l/2) 

-  ^  A 

+  S5-(Q»,0i+.,&_,rt)-A5-((3i+„i3t+„a+,/,)]'  +  AT  (11) 

In  Eq.  (11),  is  an  approximation  to  Whereas  the  superscript  n  represents  the 

current  time  level,  the  superscript  p  represents  the  iteration  number  for  the  current  time 
level.  When  p=0,  then  =  Q",  and  when  Eq.  (11)  is  iterated  to  convergence  at  a  given 
time  step,  then  In  the  present  study,  in  which  steady  state  solutions  were 


sought,  two  iterations  per  time  step  were  executed  in  all  cases.  As  discussed  by  Rai  and 
Chakravarthy  (1986),  one  primary  advantage  of  the  iterative  approach  is  that  the  stability 
of  the  system  of  equations  is  enhanced,  allowing  larger  time  steps  to  be  taken  than  with  non¬ 
iterative  approaches.  One  disadvantage  of  the  iterative  approach  is  the  increased  computer 
memory  requirement  of  needing  to  store  both  and  Q^. 

4.4  Boundary  and  Initial  Conditions 

Non-Zonal  Boundaries 

As  already  mentioned,  the  inflow  boundary  conditions  for  Zone  1  are  generated  by  a 
PNS  computation  which  provides  the  flow  conditions  at  an  axial  position  upstream  from 
the  projectile  base.  The  outer  boundary  points  of  Zones  1  and  2  are  assigned  free-stream 
values  and  located  beyond  the  bow  shock,  which  is  captured  within  the  interior  of  the 
computational  domain.  The  downstream  outflow  boundary  (Zone  2)  is  specified  using  a  first 
order  extrapolation  of  the  dependent  flow  variables.  Circumferential  boundary  conditions 
are  applied  in  both  zones  by  employing  symmetry  principles  on  the  two  planes  corresponding 
to  roll  angles  4>  =  0°  and  ^  =  30° . 

Adiabatic  boundary  conditions  are  applied  along  all  solid  surface  boundaries.  The  no¬ 
slip  boundary  condition  is  enforced,  and  a  zero  normal  pressure  derivative  is  imposed.  First 
order  implicit  boundary  conditions  are  implemented  (Rai  1989). 

Zone  2  inner  boundary  points  are  assigned  values  that  depend  upon  whether  a  miniature 
sting  or  an  axis  boundary  is  prescribed  within  the  flow  domain.  The  miniature  sting  is 
treated  as  a  solid,  adiabatic  surface.  The  axis  boundary  is  handled  using  the  following  first 
order  boundary  conditions  which  are  applicable  to  zero-yaw  flight  conditions: 

Pi  =  P2  U\  =  U-i  Ui  =  0  =  0  Pi  =  P2  (12) 

in  which  the  subscripts  1  and  2  refer  to  the  axis  boundary  point  and  first  interior  grid  point, 
respectively,  along  a  radial  grid  line.  For  the  three-dimensional  flow  situation  in  this  study, 
the  above  relationships  yield  multi-valued  axis  grid  points.  The  mean  values  of  each  of  these 
variables  are  calculated  and  used  at  all  coincident  axis  points;  the  fifth  dependent  variable, 
e,  is  then  determined  from  the  perfect  gas  relationship. 

Zonal  Boundaries 

In  the  present  grid  configuration,  the  zonal  boundary  interface,  consisting  of  two  over¬ 
lapping  boundaries,  is  located  at  the  projectile  base.  Specifically,  the  two  zonal  boundaries 
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correspond  to  the  Zone  1  outflow  plane  and  the  Zone  2  inflow  plane.  Both  zonal  boundaries 
are  perpendicular  to  the  projectile  axis  and  overlap  each  other  by  one  or  more  grid  points, 
as  shown  in  Figures  3  and  4.  The  Zone  2  zonal  boundary  is  coplanar  with  the  base,  the 
fin  trailing  edge,  and  the  next-to-last  downstream  grid  plane  of  Zone  1.  The  Zone  1  zonal 
boundary  is  located  downstream  from  the  base  and  fin  trailing  edge  and  is  constructed  to 
be  coplanar  with  either  the  first  or  one  of  the  first  upstream  interior  grid  planes  of  Zone  2. 
Since  the  zonal  boundary  points  do  not  match  those  of  the  overlapping  grid,  it  is  necessary 
to  interpolate  the  values  of  the  dependent  variables  within  each  zonal  boundary  plane. 

In  the  present  study,  an  explicit,  triangular,  Cartesian  interpolation  procedure  is  used. 
The  procedure  is  nonconservative  and  depends  solely  on  the  grid  geometries,  and  so  the 
interpolation  coefficients  associated  with  each  zonal  boundary  point  are  only  calculated  at 
the  onset  of  the  computation.  The  procedure  is  illustrated  in  Figure  5,  in  which  P  is  a  zonal 
boundary  point  whose  dependent  variables  are  to  be  defined  using  points  1,  2,  and  3  from  the 
neighboring  zone.  Point  P  is  assumed  to  be  coplanar  with  and  encompassed  by  points  1,  2, 
and  3.  In  the  current  approach,  one  of  the  three  points  is  point  P’s  closest  neighbor.  Point  O 
is  defined  by  points  1,  2,  and  3  so  that  it  forms  the  origin  of  a  right-handed,  orthogonal,  sub¬ 
coordinate  system  defined  by  unit  vectors  ei,  62,  and  63.  Let  /  represent  a  scalar  quantity 
{p  or  e,  for  example);  then  its  value  at  point  P  is  denoted  fp  and  is  found  from 

fp  =  Ci/i  -f-  €2/2  +  Ca/a  (13) 

In  this  equation,  /i,  /a,  and  /a  are  the  current,  known  values  of  /  at  points  1,  2,  and  3, 
respectively;  and  Ci,  C2,  and  Ca  are  interpolation  coefficients  given  by 

C„  =  fp  •  e„  for  n  =  l,2,  and  3  (14) 

in  which  the  vector  fp  represents  the  vector  extending  from  point  0  to  point  P. 

The  interpolation  procedure  is  not  applied  to  those  Zone  2  boundary  grid  points  which 
axe  surface  points,  i.e.,  that  constitute  the  projectile  base  or  fin  trailing  edge.  Those  points 
are  instead  flagged  at  the  onset  of  the  computation  and  are  updated  using  the  adiabatic, 
no-slip  boundary  procedure  already  mentioned. 

Initial  Conditions 

Initial  conditions  had  to  be  generated  anew  at  least  once  at  each  Mach  number.  For  a 
time-iterative  problem  of  this  magnitude,  the  initial  conditions  can  significantly  impact  the 
computer  time  needed  for  the  solution  to  approach  a  steady  state.  Zone  1  initial  conditions 
were  straightforward  because  of  the  availability  the  P  NS  solutions.  Zone  2  initial  conditions 
proved  to  be  more  sensitive,  since  free  stream  or  outer  flow  initial  values  caused  major 
numerical  instabilities  to  arise  in  almost  all  cases  when  they  were  used. 
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An  alternate  initialization  procedure  was  implemented  as  follows:  The  Zone  1  dependent 
variables  at  all  axial  locations  are  assigned  values  from  the  inflow  plane  using 


—  flj,k  (15) 

The  Zone  2  inflow  dependent  variables  are  obtained  by  applying  the  zonal  boundary 
interpolation  scheme  at  all  zonal  boundary  points,  as  previously  described.  The  cylinder 
base  and  fin  trailing  edge  surface  grid  points  are  assigned  zero  velocity,  free  stream  total 
energy,  and  density  values  from  the  perfect  gas  relation.  Next,  the  Zone  2  outflow  variables 
are  assigned  free  stream  values.  Finally,  the  interior  grid  points  are  assigned  values  by 
interpolating  between  the  Zone  2  inflow  and  outflow  grid  planes. 


5.  COMPUTATIONAL  DATA  SETS  AND  GRIDS 
5.1  Overview 

Results  from  six  CFD  solutions  are  discussed  herein  and  are  introduced  in  Table  1.  All 
but  one  of  the  data  sets  were  generated  on  a  Cray  X-MP/48  computer  (formerly  located  at 
Aberdeen  Proving  Ground  [APG]),  using  one  processor  and  limited  to  2  million  words  (2  Mw) 
of  memory.  Upon  transfer  from  NASA  Ames  Research  Center,  Moffett  Field,  California,  the 
UWIN  code  was  configured  to  load  into  core  memory  and  process  one  and  only  one  zone  at 
a  time  to  increase  the  number  of  allowable  grid  points  in  each  of  the  two  zones.  The  most 
recent,  and  final,  solution  data  set  of  the  study  (Solution  3-B)  was  generated  on  a  Cray-2 
computer  located  at  APG,  using  approximately  11  Mw  of  the  256  Mw  of  available  memory. 
The  data  sets  generated  on  the  Cray  X-MP  used  about  200  CPU  hours  each,  and  the  data 
set  generated  on  the  Cray-2  used  about  1200  CPU  hours. 

Tables  2  and  3  list  the  values  of  the  important  grid  parameters  that  were  used  for  each 
solution,  and  these  parameters  are  illustrated  in  Figure  4.  Examples  of  actual  grids  are 
shown  in  Figures  6,  7,  and  8.  Grid  clustering  is  evident  in  several  areas  and  in  all  three 
coordinate  directions.  In  any  given  direction,  the  magnitude  of  the  grid  point  clustering 
factor,  Fi,  was  kept  within  ±  20%  (but  usually  within  ±  15%)  for  all  solutions  presented 
here.  The  clustering  factor  is  defined  here  as 

F,  =  X  100%  (16) 

in  which  As/  and  As/+i  are  adjacent  grid  point  spacings  between  grid  indices  /  and  /  -f  1  in 
any  of  the  three  computational  coordinate  directions. 
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Table  1.  Base  Flow  Computational  Solution  Data  Sets 


Data  Set 

Moo 

Computer 

Data  Set 
Size  (Mb) 

Job 

Size  (Mb) 

CPU  Seconds 
per  Time  Step 

3-A 

3.0 

Cray  X-MP 

10.5 

1.8 

21. 

3-B 

3.0 

Cray  2 

66.0 

11.3 

147. 

4-A 

4.0 

Cray  X-MP 

3.6 

0.6 

6. 

4-B 

10.5 

1.8 

18. 

4-C 

4.0 

Cray  X-MP 

10.5 

1.8 

21. 

5-A 

5.0 

Cray  X-MP 

10.5 

1.8 

21. 

Table  2.  Grid  Parameters  -  Zone  1 


(All  distances  and  lengths  in  calibers) 


Data  Set 

Grid  Points 

axial  xcircx  out 

^in 

fin 

^^wall 

xlO® 

3-A 

10  X  35  X  45 

13.36 

.064 

5.8 

3-B 

30  X  35  X  45 

13.36 

.00036 

5.8 

8  X  25  X  45 

13.36 

.082 

72.0 

4-B 

10  X  35  X  45 

13.36 

.064 

5.8 

4-C 

10  X  35  X  45 

13.36 

.064 

5.8 

5-A 

10  X  35  X  45 

13.36 

.064 

Table  3,  Grid  Parameters  -  Zone  2 


(All  distances  and  lengths  in  calibers) 


Grid 

axial  xcircx  out 

^out 

^^base 

Sting 

Radius 

Aro 

Ari 

Arj 

No.  T.E. 

Points 

mm 

30  X  35  X  75 

17.91 

.0144 

none 

.0043 

.0043 

.022 

10 

54  X  71  X  160 

17.91 

.00036 

none 

.0043 

.0043 

.072 

15 

4-A 

24  X  25  X  45 

18.77 

.082 

none 

none 

2  to  4 

4-B 

25  X  35  X  90 

16.10 

.0144 

.00072 

.0036 

.0022 

.022 

10 

4-C 

30  X  35  X  75 

17.91 

.0144 

none 

.022 

10 

30  X  35  X  75 

17.91 

.0144 

none 

.0043 

.0043 

.022 

10 

11 


Figure  6  shows  a  (/>  =  0°  grid  slice  of  the  PNS  grid  (in  which  only  a  few  of  the  ax¬ 
ial  stations  in  the  region  are  shown),  as  well  as  the  Zone  1  and  Zone  2  grids  of  Solution 
4-C.  Without  exception,  the  gridding  is  comprised  of  2-D  computational  grid  planes  aligned 
perpendicular  to  the  projectile  axis.  In  Zone  1,  axial  grid  point  spacing  is  constant  except 
for  Solution  3-B,  where  points  were  clustered  near  the  base.  In  Zone  2,  axial  grid  point 
clustering  is  prescribed  near  the  projectile  base  and  radial  grid  point  clustering  is  prescribed 
at  three  radial  locations:  (1)  the  axis  (or  miniature  sting),  (2)  the  cylinder  edge,  and  (3)  the 
fin  tip. 

Figure  7  shows  a  cross  section  of  the  Zone  1  grid  of  Solution  3-B.  Each  Zone  1  grid 
cross  section  is  a  planar  elliptic  grid  constructed  using  the  same  approach  as  the  PNS  grid. 
Grid  points  are  equally  distributed  in  the  axial  direction,  with  the  exception  of  Solution 
3-B,  which  clusters  the  points  axially  near  the  base.  Grid  points  are  also  radially  clustered 
near  the  forebody  and  fin  surfaces  to  adequately  resolve  the  attached  boundary  layer.  The 
distance  from  the  wall  to  the  first  interior  grid  point  in  Zone  1  was  determined  from  the  PNS 
solutions,  which  used  the  adaptive  grid  procedure  discussed  by  Sturek  and  Schiff  (1981).  The 
value  of  the  boundary  layer  coordinate,  y'^ ,  at  the  first  point  from  the  surface  is  maintained 
in  the  range  from  2  to  5. 

Figure  8  shows  a  cross  section  of  the  Zone  2  grid  of  Solution  3-B.  Each  Zone  2  grid  cross 
section  is  algebraically  constructed.  The  Zone  2  grid  is  designed  so  that  one  grid  line  on  the 
inflow  boundary  extends  piecewise  continuously  from  the  axis  to  the  outer  boundary  and  is 
coincident  with  the  fin  planform  edge,  at  least  along  the  straight  section.  Such  construction 
aligns  the  grid  with  most  of  the  fin  planform  edge,  allowing  the  fin  boundary  layer  to  be 
better  resolved  as  it  separates  from  the  trailing  edge  and  forms  a  free  shear  layer.  Grid 
points  are  circumferentially  clustered  between  this  fin  edge  and  (/>  =  30°.  Grid  points  are 
equally  distributed  across  the  thickness  of  the  trailing  edge,  with  the  exception  of  Solution 
3-B,  which  clusters  the  trailing  edge  points  near  the  fin  planform  edge. 

5.2  Details  of  Data  Sets 

The  discussion  of  the  individual  CFD  data  sets  begins  at  Mach  4,  the  free  stream  velocity 
of  the  initial  problem  setup  upon  transfer  from  NASA  Ames.  The  data  sets  are  described 
here  in  mostly  chronological  order  by  first  discussing  the  Mach  4  cases,  followed  by  the  Mach 
3  cases,  then  the  Mach  5  Ccise. 
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Solution  4- A.  at  Mach  4,  is  the  first  solution  generated  in  the  study  that  includes  the 
projectile  base  region.  It  is  chronologically  the  first  solution  to  be  generated  of  the  six 
solutions  presented  here.  The  grid  resolution  for  Solution  4-A,  as  apparent  from  Table  3,  is 
considerably  less  than  the  other  five  CFD  solutions.  The  Zone  2  outvirard  grid  lines  of  Solution 
4-A  are  oriented  along  constant  <l>  values,  as  opposed  to  being  aligned  with  the  fin  planform 
edge.  Solution  4-A  also  includes  a  miniature  sting,  since  the  axis  boundary  condition  had 
yet  to  be  implemented.  Solution  4-A  was  successfully  initialized  using  free-stream  values  for 
the  interior  grid  points. 

Solution  4-B.  also  at  Mach  4,  is  the  first  of  the  study  to  radially  cluster  points  at  three 
radial  locations:  the  axis,  the  projectile  cylinder  edge,  and  the  fin  tip.  It  is  also  the  first 
solution  to  circumferentially  align  the  Zone  2  grid  with  the  fin  planform  edge,  as  was  illus¬ 
trated  in  Figure  6.  Like  Solution  4-A,  it  includes  a  miniature  sting.  The  outflow  boundary  of 
Solution  4-B  was  moved  upstream  by  2.67  cal  compared  to  that  of  Solution  4-A  in  an  effort 
to  save  grid  points.  Solution  4-B  was  initialized  using  the  alternate  procedure  described 
earlier. 

Solution  4-C.  the  third  and  final  solution  at  Mach  4,  is  the  first  solution  of  the  study 
to  include  the  axis  boundary  condition.  The  Zone  2  radial  grid  spacing  and  distribution 
were  modified  slightly  compared  to  Solution  4-B.  Also,  the  outflow  boundary  was  moved 
downstream  by  1.81  cal  compared  to  that  of  Solution  4-B  since  it  was  noticed  that  Solution 
4-B  had  regions  of  subsonic  velocity  in  the  outflow  plane.  Solution  4-C  was  initialized  by 
interpolating  the  dependent  variables  from  Solution  4-B. 

Solution  3-A.  at  Mach  3,  includes  all  the  features  of  Solution  4-C,  with  the  exception  of 
minor  differences  in  Zone  2  grid  spacing  and  distribution  (see  Table  3).  Solution  3-A  was 
initialized  using  the  alternate  procedure  described  earlier. 

Solution  3-B.  also  at  Mach  3,  is  chronologically  the  last  solution  to  be  generated  during 
the  study.  It  has  almost  six  times  more  grid  points  in  Zone  2  (over  600,000  total)  than  the 
other  data  sets.  The  additional  grid  points  were  used  to  increase  the  near  wall  resolution 
in  the  axial  and  circumferential  directions.  The  number  of  grid  points  across  the  trailing 
edge  half-thickness  was  15,  and  Solution  3-B  was  the  only  data  set  to  circumferentially 
cluster  these  trailing  edge  points  along  the  fin  planform  edge.  Solution  3-B  was  initialized 
by  interpolating  the  dependent  variables  from  Solution  3-A. 

Solution  5-A.  the  only  solution  at  Mach  5,  has  the  same  modeling  and  grid  features  as 
Solutions  3-A  and  4-C.  Solution  5-A  was  initialized  using  the  alternate  procedure  described 
earlier. 
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6.  RESULTS 


6.1  Drag  Comparison.  The  range- acquired  drag  coefficients  (see  Table  4)  were  ob¬ 
tained  from  unpublished  firings  conducted  at  the  indoor  ARL  transonic  range  facility  at 
APG  during  1979.  The  14  shots  listed  here  were  fired  with  an  empty  tracer  cavity.  The 
zero-yaw  drag  coefficients  were  produced  from  standard  yaw-drag  fits  of  the  equations  of 
motion  from  Murphy  (1963).  The  total  drag  coefficient,  Cd-,  of  a  projectile  at  some  instant 
in  time  is  defined  as 

Cd  =  Cdo  +  Cofi^^  (17) 

in  which  Cdq  is  the  zero-yaw  drag  coefficient,  i®  ^rag  coefficient,  and  6  is  the 

angle  of  attack.  The  set  of  14  fitted  values  of  Cd  and  6  were  used  to  make  a  least  squares 
fit  to  the  equation 

Cd  =  Ai  +  +  (18) 

in  which  values  of  Ai  and  Aj  (coefficients  of  no  interest  here),  as  well  as  Cd  were  obtained. 
The  Mach  number  for  each  shot  is  Moo-  The  resulting  value  of  14.6  for  Cd^^  was  then  used 
to  repeat  the  yaw-drag  fits  to  obtain  Cdo- 


Table  4.  Range-Acquired  Zero- Yaw  Drag  Coefficients 


Shot  No. 

Moo 

Cdo 

Shot  No. 

Afoo 

Cdo 

16405 

4.295 

0.295 

16417 

4.113 

0.305 

16423 

4.289 

0.293 

16401 

3.873 

0.324 

16404 

4.281 

0.302 

16402 

3.869 

0.327 

16424 

4.262 

0.300 

16403 

3.759 

0.342 

16421 

4.237 

0.303 

16419 

3.447 

0.372 

16422 

4.223 

0.305 

16418 

3.437 

0.368 

16416 

4.123 

0.306 

16420 

3.015 

0.420 

Table  5  lists  the  computed  zero-yaw  drag  coefficients  and  various  components  for  the 
six  CFD  data  sets.  As  mentioned  earlier,  the  computed  zero- yaw  drag  coefficient,  Cdo->  is 
defined  as 

Cdo  =  Cdf  +  Cdb  (19) 

in  which  Cdp  is  the  forebody  drag  contribution  (which  includes  the  forebody  and  fin  con¬ 
tribution,  i.e.,  the  attached  flow)  obtained  from  the  PNS  code,  and  Cdq  is  the  base  drag 
contribution  (which  includes  the  cylinder  base  and  fin  trailing  edge  drag  contribution)  ob¬ 
tained  from  the  UNS  code.  The  computed  drag  coefficients  are  calculated  by  using  the  force 
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Table  5.  Computed  Zero- Yaw  Drag  Coefficients 


Moo 

Data  Set 
(UNS) 

Cdb 

(UNS) 

Cdo 

{PB  !  Poo)  avg 

3.0 

3-A 

0.123 

0.363 

0.463 

99 

3-B 

0.376 

4.0 

0.222 

4-A 

0.047 

0.269 

0.636 

r 

99 

4-B 

0.084 

0.306 

0.349 

99 

99 

4-C 

0.082 

0.304 

0.364 

5.0 

0.202 

5-A 

0.047 

0.249 

0.433 

acting  on  the  entire  base  and  six  fins  of  the  projectile,  consistent  with  the  range-acquired 
drag  coefficients.  At  Mach  4,  the  computed  base  drag  of  Solution  4-A  (the  coarsest  grid  of 
the  six  CFD  data  sets)  is  about  56%  of  that  of  Solutions  4-B  or  4-C.  The  computed  base 
drag  of  Solutions  4-B  and  4-C  themselves  differs  by  less  than  3%.  At  Mach  3,  the  base  drag 
of  Solution  3-B  (the  finest  grid  of  the  six  CFD  data  sets)  is  about  15%  greater  than  that  of 
Solution  3- A. 

Figure  9  shows  the  comparison  of  computed  and  range-acquired  zero-yaw  drag  coeffi¬ 
cients.  The  results  of  Solutions  3-B,  4-C,  and  5-A  are  combined  with  the  PNS  results  (also 
shown)  in  the  figure.  The  computed  drag  underpredicts  the  range-acquired  drag  by  about 
9%  and  5%,  at  Mach  3  and  4,  respectively.  Since  no  firings  were  conducted  above  Mach  4.3, 
no  comparison  can  be  made  for  the  Mach  5  case.  The  computed  base  drag  comprises  about 
33%,  26%,  and  19%  of  the  computed  total  drag  at  Mach  3,  4,  and  5,  respectively. 

Though  the  drag  comparison  serves  as  an  initial  indicator  of  the  accuracy  of  the  CFD 
approach,  more  precise  validation  of  the  predicted  flow  field  is  needed.  The  focus  of  future 
research  in  this  area  should  be  the  validation  of  base  and  fin  trailing  edge  pressure  predictions 
against  wind  tunnel  measurements  such  as  those  reported  by  Berner  et  al.  (1992)  or  cited 
by  Moore  et  al.  (1993).  While  no  further  quantitative  comparisons  are  made  in  this  report, 
some  details  of  the  computed  base  flow  structure  are  presented. 


6.2  Base  Pressure.  Table  5  also  includes  the  average  base  pressure,  (ps/poo)ova,  for 
each  of  the  six  CFD  solutions.  The  average  base  pressure  consists  of  the  cylinder  base  and 
fin  trailing  edge  contributions.  A  comparison  of  Solutions  3-A,  4-B,  and  5-A,  which  basically 
have  the  same  grid  resolution,  reveals  that  the  computed  average  base  pressure  of  Mach  5 
is  less  than  that  of  Mach  3  but  greater  than  that  of  Mach  4.  Such  non-monotonic  behavior 
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with  respect  to  free  stream  Mach  number  is  unexpected,  at  least  for  axisymmetric  base  flows 
as  well  as  2-D  trailing  edge  flows.  The  result  does  not  specifically  eliminate  the  possibility 
that  interaction  between  the  axisymmetric-like  cylinder  base  and  2-D-like  trailing  edge  could 
produce  such  behavior  but  should  be  viewed  as  a  possible  inconsistency  between  solutions. 
The  result  may  be  attributable  to  inadequate  grid  resolution,  modeling  aissumptions,  or 
insufficient  convergence  of  one  or  more  of  the  CFD  solutions.  A  conclusive  interpretation 
cannot  be  made  without  quantitative  pressure  measurements  with  which  to  compare. 

Figures  10  through  14  show  the  computed  base  pressure  distributions  along  roll  angles 
<f>  =  0°  and  (j)  =  30°  for  each  of  the  six  CFD  solutions.  At  Mach  4,  Solution  4-A  has 
much  higher  base  pressures  than  Solutions  4-B  and  4-C,  almost  certainly  attributable  to  its 
comparatively  sparse  grid.  Solution  4-B  has  a  lower  base  pressure  on  the  cylinder  portion 
than  Solution  4-C,  but  the  trend  is  reversed  on  the  fin  trailing  edge,  and  the  integrated  effect 
leads  to  the  similar  base  drag  values  for  the  two  solutions,  as  already  shown.  At  Mach  3, 
this  same  pressure  reversal  on  cylinder  base  and  fin  trailing  edge  is  found  between  Solutions 
3- A  and  S-B. 

All  the  CFD  solutions  except  Solution  5-A  show  higher  pressure  along  the  fin  trailing 
edge  than  on  the  cylinder  base.  Along  (f>  =  30°,  all  six  CFD  solutions  show  at  least  a  slight 
undershoot  in  pressure  near  the  cylinder  edge,  indicative  of  a  possible  need  for  improved  grid 
resolution  there.  The  same  observation  may  apply  along  <^  =  0°  near  the  fin  tip.  For  each 
CFD  solution,  small  differences  between  pressure  distributions  along  <^  =  0°  and  (/>  =  30°  are 
apparent,  but  the  details  are  inconsistent;  in  some  cases,  higher  pressures  exist  along  <^  =  0° 
than  ^  =  30°,  and  in  other  cases,  the  opposite  is  true.  Of  all  the  pressure  distributions 
shown  here,  that  of  Solution  3-B  has  the  most  apparent  saddle  points,  probably  attributable 
to  its  finer  grid  resolution  compared  to  the  other  five  CFD  solutions. 

6.3  Flow  Structure 

Mach  4  /  Mach  4.3 

Figure  15  is  a  reproduced  shadowgraph  mosaic  of  an  M735  projectile  at  Mach  4.3,  taken 
at  the  ARL  transonic  range  facility  during  a  nosetip  boundary  layer  transition  study  reported 
by  Sturek,  Kayser,  and  Mylin  (1983).  The  shadowgraph,  discovered  late  in  the  present  study, 
is  included  here  to  at  least  qualitatively  compare  with  Solution  4-C,  which  is  at  Mach  4.0. 
The  film  surface  is  estimated  to  be  parallel  to  roll  angle  (j)  =  24°  within  ±3°. 

Two  main  features  apparent  in  the  base  region  of  the  shadowgraph  are  (1)  a  turbulent 
flow  region  extending  practically  from  tip-to-tip  of  opposite  fins  and  (2)  two  or  three  shocks 
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emanating  from  an  axial  location  near  the  cylinder  base  above  and  below  the  projectile  axis. 
The  shocks  are  oriented  at  approximately  10°  relative  to  the  axis.  The  two  shocks  nearest 
to  the  axis  almost  appear  to  intersect  the  axis  at  about  1  cal  down  stream  from  the  base. 
From  the  shadowgraph  alone,  an  observer  cannot  conclusively  determine  whether  all  the 
shock  images  in  the  cylinder  base  region  represent  structures  that  are  coplanar  with  the 
<j>  =  24°  ±  3°  plane,  as  opposed  to  being  parallel  to  it. 

Figures  16  and  17  show  computed  Mach  contours  for  Solution  4-C  in  the  ^  =  0°  and 
(f>  =  30°  planes.  In  the  (^  =  0°  plane,  a  well  defined  velocity  wake  region  is  apparent  from 
tip-to-tip  of  opposite  fins.  The  velocity  wake  region  is  also  well-defined  in  the  <j>  =  30° 
plane,  contracting  to  a  minimum  width  of  about  3/4  cal  at  about  3/4  cal  downstream, 
and  expanding  to  a  maximum  width  of  about  1.5  cal  at  about  3  cal  downstream.  Both 
roll  angles  show  a  small  core  (widths  on  the  order  of  1/10  cal)  of  low  Mach  number  flow 
extending  downstream  along  the  axis.  Both  roll  angles  also  show  the  width  of  this  core  to 
noticeably  increase  about  3  cal  downstream  from  the  base,  especially  in  the  ^  =  0°  plane. 

Figures  18  and  19  show  computed  pressure  contours  for  Solution  4-C  in  the  ^  =  0° 
and  <t>  =  30°  planes.  Some  qualitative  similarity  in  the  location  and  orientation  of  shocks 
can  be  discerned  between  Solution  4-C  and  the  shadowgraph.  In  the  CFD  solution,  the 
shock  structures  more  closely  resemble  recompressions  since  they  are  highly  smeared  (and 
somewhat  step-like  in  shape),  probably  attributable  to  inadequate  grid  resolution.  As  in 
the  shadowgraph,  the  computed  shock  structures  are  oriented  at  approximately  10°  relative 
to  the  axis.  They  are  more  pronounced  in  the  ((>  =  30°  plane  than  in  the  =  0°  plane, 
and  those  nearest  the  axis  extend  no  farther  upstream  than  about  1  cal  from  the  base.  In 
addition,  the  pressure  contours  in  both  planes  show  a  region  of  relatively  compressed  flow 
along  the  axis  beginning  0.5  cal  downstream  of  the  base  and  extending  about  1.5  cal  farther 
downstream. 

In  the  previous  figures,  no  major  discontinuities  are  apparent  in  the  contour  lines  that 
intersect  the  zonal  boundary  interface.  On  a  smaller  scale,  however,  it  is  noted  that  some 
minor  numerical  disturbances  were  introduced  into  all  from  the  flow  solutions  near  the  curved 
edges  of  the  base  and  fin  trailing  edge.  Two  likely  causes  of  these  disturbances  are  (1)  the  use 
of  the  current  zonal  boundary  interpolation  scheme  in  conjunction  with  such  relatively  sparse 
grid  resolution  in  the  vicinity  of  the  solid  boundaries,  and  (2)  the  use  of  fully  turbulent  flow 
conditions  in  Zone  1  and  fully  laminar  flow  conditions  in  Zone  2  introducing  an  additional 
flow  discontinuity  to  numerically  capture. 
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Mach  3 


Of  the  six  CFD  solutions  discussed  in  this  report,  Solution  3-B,  at  Mach  3,  has  the  most 
grid  points  and  best  spatial  resolution.  Figures  20  through  22  show  Mach  contours,  pressure 
contours,  and  velocity  vectors,  respectively,  in  the  ^  =  0°  and  (j)  =  30°  half-planes  of  Zone 
2.  In  these  figures,  the  two  half-planes  are  presented  perpendicular  to  each  other,  with  the 
^  =  0°  half-plane  above  the  axis.  These  figures  show  a  closer  view  of  the  flow  near  the  base 
than  the  Mach  4  case  just  shown. 

In  the  (/>  =  0°  half-plane.  Figures  20  and  21  show  one  shock  structure  emanating  from  a 
location  near  the  axis  at  an  angle  of  approximately  30°  relative  to  the  axis.  This  is  in  contrast 
to  the  Mach  4  case  just  shown,  which  had  two  shock  structures  in  each  half-plane.  As  in  the 
Mach  4  case,  smearing  and  step-like  behavior  is  apparent  in  the  Mach  and  pressure  contours 
that  comprise  this  computed  shock  structure.  The  shock  structure  practically  intersects  the 
axis  at  a  location  about  1  cal  downstream  from  the  base.  There,  a  compressed  low  Mach 
number  region  of  axial  length  about  1/2  cal  and  extending  1/10  cal  in  the  radial  direction 
is  apparent  in  both  half-planes.  Downstream,  the  contours  intersect  the  axis  in  a  more 
perpendicular  fashion  than  do  those  of  the  Mach  4  solution  already  shown. 

In  Figure  22,  the  flow  directly  downstream  from  the  fin  trailing  edge  moves  strongly 
toward  the  axis  and  meets  outwardly  moving  flow  from  the  cylinder  base  region.  Figures  20 
and  22  show  that  the  flow  from  these  two  regions  is  separated  by  a  weak  shear  layer  oriented 
about  30°  relative  to  the  axis. 

In  the  ^  =  0°  half-plane.  Figure  22  shows  a  small  region  of  reverse  flow  (that  is,  in  which 
the  axial  velocity  component  is  oriented  upstream)  within  1/2  cal  downstream  from  the  fin. 
Within  the  cylinder  base  region  itself,  no  reverse  flow  is  computed  to  occur  within  1/2  cal 
from  the  base,  indicating  that  the  fluid  is  entering  the  ^  =  0°  half-plane  near  the  base. 
Approximately  1  cal  farther  downstream  in  the  <^  =  0°  half-plane,  reverse  flow  is  apparent, 
indicating  that  the  fluid  is  exiting  the  ^  =  0°  half-plane. 

In  the  (j)  =  30°  half-plane.  Figure  22  shows  reverse  flow  within  1  /4  cal  from  the  base.  The 
reverse  flow  is  part  of  a  recirculation  zone  whose  center  marks  the  upstream  edge  of  a  shear 
layer.  The  shear  layer,  located  approximately  1/4  cal  from  the  axis,  separates  the  reversed 
and  non-reversed  flow  regions  and  extends  farther  downstream.  The  free  shear  layer  that 
enters  the  (f>  =  30°  plane  from  the  cylinder  surface  is  shown  to  largely  dissipate  by  the  time 
the  flow  reaches  1  cal  downstream. 

Figure  23  shows  computed  axial  velocity  contours  for  Solution  3-B  at  six  cross-sectional 
planes  downstream  from  the  base.  The  axial  velocity,  u,  is  normalized  by  the  free  stream 
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speed  of  sound,  Goo-  The  farthest  downstream  plane  is  located  approximately  1  cal  from  the 

ase.  For  purposes  of  orientation,  the  projectile  axis  is  plotted  as  a  thick  solid  line.  In  the 

figure,  black  contour  lines  represent  negative  axial  velocity  (i.e.,  reversed  flow)  The  figure 

clearly  shows  a  core  of  reversed  flow  embedded  within  the  computational  domain.  The  center 

of  the  core  is  located  on  neither  the  ^  =  0°  nor  the  <t>  =  30°  plane.  No  reverse  flow  occurs 

anywhere  m  the  computed  flow  field  downstream  from  the  most  downstream  plane  in  the 
figure. 

Figure  24  shows  computed  crossflow  velocity  vectors  for  Solution  3-B  within  a  cross 
^tional  plane  downstream  from  the  base.  This  plane  is  located  approximately  3/8  cal 
downstream  from  the  base  and  corresponds  to  the  fourth  downstream  plane  from  the  base 
shown  m  the  previous  figure.  The  plot  includes  a  shaded  area  which  provides  a  magnified 
view  of  the  reversed  flow  region  from  the  previous  figure.  At  this  axial  location,  it  is  apparent 
tha  the  crossflow  direction  within  the  reverse  flow  region  is  primarily  outward  from  the  axis, 
While  the  crossflow  direction  elsewhere  is  primarily  inward  toward  the  axis. 


7.  CONCLUSION 


CFD  study  of  finned  KE  projectile  base  flow  has  been  presented  for  supersonic  flight 
conditions.  A  parabolized  Navier-Stokes  technique  was  used  to  compute  the  attached  flow 
over  the  forebody  and  fins;  an  unsteady  Navier-Stokes  technique  was  used  to  compute  the 
flow  in  the  base  region.  Several  CFD  solutions  for  free  stream  Mach  numbers  3,  4,  and  5 
were  presented,  and  the  modeling  approach  and  computational  requirements  were  discussed 

The  computed  zero-yaw  drag  was  compared  to  range-acquired  values,  and  the  flow  structure 
was  bneny  examined. 

r  ^  Y  »”  “>'=  order  of  hundreds  of  hours  on  a 

Cr  j  X-MP  computer.  The  largest  solution,  which  had  a  six-fold  increase  in  the  number  of 

^.d  points  (more  than  600,000)  compared  to  the  others,  used  more  than  1000  hours  on  a 
Cray-2  computer.  All  of  the  solution  data  sets  are  raluable,  not  only  because  they  represent  a 

Tfut  rFr‘T"‘ 

lor  future  CFD  work  in  this  area. 

While  the  present  two-zone  approach  provides  a  good  starting  point  for  the  problem,  a 
mu  tiple  zone  strategy  would  provide  more  flexibility,  especially  for  more  complex  config- 

Lnlr  r  As  discussed  within,  several 

aspects  of  the  CFD  results  strongly  suggest  that  a  further  increase  in  grid  resolution  may  be 
required  to  resolve  the  flow  structures  in  certain  regions  of  the  flow  field. 
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The  computed  zero-yaw  drag  of  the  projectile  was  shown  to  underpredict  the  range- 
acquired  drag  by  about  9%  at  Mach  3  and  5%  at  Mach  4.  The  many  modeling  simplifica¬ 
tions  used  in  the  current  CFD  approach  which  could  contribute  to  differences  between  the 
computed  and  range-acquired  drag  were  discussed.  Though  the  drag  comparison  serves  as 
an  initial  indicator  of  the  accuracy  of  the  CFD  approach,  it  is  concluded  that  precise  vali¬ 
dation  of  finned  projectile  base  flow  CFD  modeling  should  focus  on  wind  tunnel  base  and 
fin  trailing  edge  pressure  measurements.  Future  work  in  this  area  should  also  include  the 
investigation  of  finless  models. 

Overall,  the  computed  base  flow  structure  is  shown  to  exhibit  complex,  fully  3-D  features. 
No  axisymmetry  exists  within  the  cylinder  base  region,  and  the  base  flow  structure  is  quite 
different  from  that  of  axisymmetric  base  configurations  that  have  been  reported  to  date.  It 
is  conceded  that  uncertainty  remains  concerning  the  adequacy  of  the  current  grid  resolution 
and  modeling  simplifications.  Much  of  this  uncertainty  is  driven  by  the  lack  of  sufficient 
computational  power  needed  to  conduct  parametric  studies  that  address  these  issues.  The 
current  facilities,  while  having  enough  core  memory,  provide  insufficient  CPU  time  per  user 
to  allow  such  an  investigation  to  be  promptly  conducted.  The  study  presented  here  serves 
as  a  stepping-stone  for  future  CFD  modeling  that  should  eventually  lead  to  a  capability  for 
predicting,  understanding,  and  controlling  KE  projectile  base  flow. 
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Figure  2.  Roll  Angle,  <j) 


KE  PROJECTILE  BASE  ELDW  MODELING 


Figure  3.  Code  Coupling  for  Finned  Projectile  Base  Flow  Modeling,  (i> 


Figure  7.  Cross-Sectional  View  of  Grid  3-B,  Zone  1,  x/d=  13.94 


Figure  8.  Cross-Sectional  View  of  Grid  3-B,  Zone  2,  x/d=13.94 
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Mach  Number 


Figure  9.  Comparison  of  Computed  and  Range- Acquired  Zero- Yaw  Drag  Coefficients 
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Figure  10.  Computed  Base  Pressures,  Mqo  =  4,  ^  =  0° 
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Figure  12.  Computed  Base  Pressures,  Moo  =  3,  </>  =  0° 
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Figure  13.  Computed  Base  Pressures,  Moo  =  3,  (^  =  30° 
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Figure  14.  Computed  Beise  Pressures,  Moo  =  b,  (f)  =  ^°  and  6  =  30° 


Figure  16.  Mach  Contours  (0.0  to  4.2  in  0.1  inci 
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(Solution  4-C),  4) 


Figure  17.  Mach  Contours  (0.0  to  4.6  in  0.1  increments),  Moo  =  4  (Solution  4-C),  4>  =  30' 


Figure  20.  Mach  Contours  (0.0  to  4.0  in  0.2  increments),  =  3  (Solution  3-B),  <f)  =  0° 
and  0  =  30° 
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Figure  21.  p/poc  Contours  (0.15  to  1.2  in  0.05  increments),  M^o  =  3  (Solution  3-B),  (f)  =  0° 
and  (p  =  30° 
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LIST  OF  SYMBOLS 


A±,B±,C± 

a 

Cd 

Cdo 

^Db 

Cdf 

Cd,2 

d 

E,F,G 

E,F,G 

I 

i,j,k 

J 

L 

M 

N 

Pr 

P 

Q 

Q 

q 

Re 

s 

i 

u,v,w 

u,v,w 


flux-split  Jacobian  matrices  of  transformed  inviscid  flux  vectors 

speed  of  sound 

total  drag  coeflBcient 

zero-yaw  drag  coeflicient 

base  and  fin  trailing  edge  drag  coefficient 

forebody  and  fin  attached  flow  drag  coefficient 

yaw  drag  coefficient 

projectile  reference  diameter 

inviscid  flux  vectors  of  transformed  gas  dynamic  equations 

first  order  numerical  inviscid  flux  vectors  of  transformed  gas  dynamic  equations 

identity  matrix 

grid  indices  in  ?/,  (  directions 
Jacobian  of  coordinate  transformation 
unit  length 
Mach  number 

Jacobian  matrix  of  viscosity  terms  in  ^  direction 

Prandtl  number  for  laminar  flow 

pressure 

vector  of  dependent  variables  of  gas  dynamic  equations 

transformed  vector  of  dependent  variables  of  gas  dynamic  equations 

dynamic  pressure 

Reynolds  number,  poo^ooL! Poo 

transformed  vector  viscosity  terms  in  ^  direction 

time 

velocity  components  in  ^,7;,  and  ^  directions 
velocity  components  in  x,  j/,  and  z  directions 


Greek  Symbols 


7 

e 

P 

P 

4> 

r 


ratio  of  specific  heats 

total  energy  per  unit  volume  of  fluid 

coefficient  of  molecular  viscosity 

density 

roll  angle 

transformed  time 

transformed  coordinates 
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Subscripts 

B 

8 

wall 

oo 


projectile  base  condition 
total  angle  of  attack 
wall  condition 
free-stream  condition 


Superscripts 

T  vector  transpose 


Note: 

Aerodynamic  force  coefficients  are  defined  as 


Force 

~qoo(v<P/4) 
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